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Abstract 


This  report  describes  the  prediction  of  appendage  and  viscous  forces  in  DRDC  At¬ 
lantic’s  ShipMo3D  library  for  ship  motions  in  waves.  Previously,  the  ShipMo3D 
library  considered  only  hull  forces  due  to  potential  flow.  Inclusion  of  appendage 
and  viscous  forces  is  essential  for  accurate  prediction  of  ship  sway,  roll,  and  yaw 
motions.  Appendage  forces  are  caused  by  added  mass,  lift,  and  viscous  effects.  Hull 
viscous  forces  are  also  important  and  are  included  in  the  present  treatment.  This  re¬ 
port  introduces  a  method  for  including  appendages  in  hull  radiation  and  diffraction 
computations  by  representing  appendages  using  dipole  panels.  Computations  for 
a  naval  destroyer  have  indicated  that  inclusion  of  appendages  in  hull  radiation  and 
diffraction  computations  has  a  negligible  effect  on  predicted  motions;  thus,  hull  and 
appendage  forces  can  be  evaluated  separately  without  affecting  accuracy  of  motion 
predictions.  Comparisons  of  predictions  with  experiments  for  a  steered  warship 
model  in  regular  waves  show  good  agreement.  Future  work  will  incorporate  hull 
lift  forces  in  ship  motion  predictions. 


Resume 


Le  present  rapport  decrit  la  prediction  des  forces  exercees  sur  les  appendices  et 
des  forces  de  viscosite  dans  la  bibliotheque  ShipMo3D  des  mouvements  des  na- 
vires  dans  les  vagues  de  RDDC  Atlantique.  Auparavant,  la  bibliotheque  ShipMo3D 
tenait  compte  uniquement  des  forces  exercees  sur  la  coque  par  l’ecoulement  po- 
tentiel.  L’ inclusion  des  forces  exercees  sur  les  appendices  et  des  forces  de  vis¬ 
cosite  est  essentielle  pour  predire  avec  precision  les  mouvements  de  deplacement 
lateral,  de  roulis  et  d’embardee  des  navires.  Les  forces  exercees  sur  les  appendices 
sont  causees  par  la  masse  ajoutee,  la  portance  et  les  effets  de  viscosite.  Les  forces 
de  viscosite  de  la  coque  sont  aussi  importantes  et  sont  indues  dans  les  presentes 
operations.  Le  present  rapport  traite  d’une  methode  pour  inclure  les  appendices 
dans  les  calculs  de  radiation  et  de  diffraction  de  la  coque  en  representant  les  appen¬ 
dices  au  moyen  de  panneaux  de  doublets.  Les  calculs  effectues  pour  un  destroyer 
indiquent  que  T inclusion  des  appendices  dans  les  calculs  de  radiation  et  de  diffrac¬ 
tion  de  la  coque  a  tres  peu  d’effet  sur  les  mouvements  prevus ;  les  forces  exercees 
sur  la  coque  et  les  appendices  peuvent  done  etre  evaluees  separement  sans  que  cela 
influe  sur  la  precision  des  predictions  de  mouvements.  Les  comparaisons  entre  les 
predictions  et  les  experiences  avec  un  modele  de  navire  de  guerre  sous  manoeuvre 
dans  des  vagues  regulieres  donnent  une  bonne  concordance.  Des  travaux  futurs  in- 
cluront  les  forces  de  portance  de  la  coque  dans  les  predictions  des  mouvements  des 
navires. 
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Executive  summary 

Introduction 

DRDC  Atlantic  is  developing  a  new  object-oriented  library  for  simulation  of  ship 
motions  in  waves.  It  is  widely  known  that  appendage  and  viscous  forces  signif¬ 
icantly  influence  ship  motions  in  sway,  roll,  and  yaw.  This  report  describes  the 
implementation  of  appendage  and  viscous  forces  in  the  ShipMo3D  library. 


Principal  Results 

Appendage  and  viscous  forces  have  been  implemented  in  the  ShipMo3D  library 
for  computations  in  both  the  frequency  domain  and  time  domain.  This  report  intro¬ 
duces  a  method  for  including  appendages  in  hull  radiation  and  diffraction  compu¬ 
tations  by  representing  appendages  using  dipole  panels.  Computations  for  a  naval 
destroyer  have  indicated  that  inclusion  of  appendages  in  hull  radiation  and  diffrac¬ 
tion  computations  has  a  negligible  effect  on  predicted  motions;  thus,  hull  and  ap¬ 
pendage  forces  can  be  evaluated  separately  without  affecting  accuracy  of  motion 
predictions.  Comparisons  of  predictions  with  experiments  for  a  steered  warship 
model  in  regular  waves  show  good  agreement.  Excellent  agreement  between  fre¬ 
quency  domain  and  time  domain  predictions  indicates  consistent  implementation 
of  appendage  and  viscous  force  terms. 


Significance  of  Results 

The  ShipMo3D  library  can  now  model  motions  in  waves  for  six  degrees  of  freedom 
for  a  ship  with  steady  speed  and  heading.  Validation  to  date  suggests  that  the  library 
can  be  used  with  confidence  for  ships  in  moderate  sea  conditions. 


Future  Plans 

The  ShipMo3D  library  is  being  extended  to  model  motions  of  a  freely  maneuvering 
ship  in  waves.  Hull  lift  forces  will  be  incorporated  in  future  ship  motion  predictions. 
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Sommaire 


Introduction 

RDDC  Atlantique  travaille  a  la  mise  au  point  d’une  nouvelle  bibliotheque  orientee  objets 
pour  la  simulation  des  mouvements  des  navires  dans  les  vagues.  II  est  reconnu  que  les  forces 
exercees  sur  les  appendices  et  les  forces  de  viscosite  influent  de  maniere  significative  sur  les 
mouvements  des  navires  dans  les  situations  de  deplacement  lateral,  de  roulis  et  d’embardee. 
Le  present  rapport  traite  de  1' implementation  de  donnees  sur  les  forces  exercees  sur  les 
appendices  et  les  forces  de  viscosite  dans  la  bibliotheque  ShipMo3D. 

Resultats  principaux 

Des  donnees  sur  les  forces  exercees  sur  les  appendices  et  les  forces  de  viscosite  ont  ete 
implementees  dans  la  bibliotheque  ShipMo3D  pour  effectuer  des  calculs  dans  le  domaine 
des  frequences  et  le  domaine  temporel.  Le  present  rapport  traite  d’une  methode  pour  inclure 
les  appendices  dans  les  calculs  de  radiation  et  de  diffraction  de  la  coque  en  representant  les 
appendices  au  moyen  de  panneaux  de  doublets.  Les  calculs  effectues  pour  un  destroyer 
indiquent  que  l’inclusion  des  appendices  dans  les  calculs  de  radiation  et  de  diffraction  de 
la  coque  a  tres  peu  d’effet  sur  les  mouvements  prevus  ;  les  forces  exercees  sur  la  coque  et 
les  appendices  peuvent  done  etre  evaluees  separement  sans  que  cela  influe  sur  la  precision 
des  predictions  de  mouvements.  Les  comparaisons  entre  les  predictions  et  les  experiences 
avec  un  modele  de  navire  de  guerre  sous  manoeuvre  dans  des  vagues  regulieres  donnent  une 
bonne  concordance. 

Importance  des  resultats 

La  bibliotheque  ShipMo3D  peut  maintenant  modeliser  des  mouvements  dans  les  vagues  sur 
six  degres  de  liberte  pour  un  navire  avec  une  vitesse  et  un  cap  constants.  Jusqu’a  maintenant, 
la  validation  laisse  supposer  que  la  bibliotheque  peut  etre  utilisee  pour  les  navires  dans  des 
conditions  de  mer  agitee. 

Travaux  ulterieurs  prevus 

La  bibliotheque  ShipMo3D  enrich ic  permet  de  modeliser  les  mouvements  d’un  navire  en 
libre  manoeuvre  dans  les  vagues.  Les  forces  de  portance  de  la  coque  seront  incorporees  dans 
des  predictions  futures  des  mouvements  des  navires. 
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1  Introduction 


DRDC  Atlantic  is  developing  software  components  for  modelling  and  simulation  of 
ships  in  waves.  Work  completed  to  date  includes  predictions  of  forces  and  motions 
in  waves  for  an  unappended  ship  hull  in  the  frequency  domain  [1]  and  time  do¬ 
main  [2],  and  modelling  of  seaways  for  prediction  of  ship  motions  [3].  This  report 
describes  the  inclusion  of  appendage  and  viscous  forces  in  the  ShipMo3D  library. 

Appendage  and  viscous  forces  must  be  considered  for  accurate  prediction  of  ship 
motions  in  sway,  roll,  and  yaw.  Schmitke  [4]  was  among  the  first  to  include  ap¬ 
pendage  and  viscous  forces  in  ship  motion  predictions.  The  accurate  prediction 
of  viscous  roll  damping  continues  to  be  a  challenge.  Himeno  [5]  provides  a  very 
good  overview  of  viscous  forces  influencing  roll  damping.  These  two  references 
provide  a  significant  foundation  for  incorporating  appendage  and  viscous  forces  in 
the  ShipMo3D  library.  The  ShipMo3D  library  introduces  a  new  capability  for  in¬ 
cluding  appendages  in  hull  radiation  and  diffraction  computations,  as  discussed  by 
McTaggart  and  Stredulinsky  [6]. 

Section  2  of  this  report  describes  the  axis  system  used  for  computing  forces  on 
a  ship  hull  and  appendages.  Section  3  describes  the  prediction  of  forces  on  ship 
appendages,  followed  by  the  prediction  of  hull  viscous  forces  in  Section  4.  The  in¬ 
clusion  of  appendages  in  hull  radiation  and  diffraction  computations  is  described  in 
Section  5.  Section  6  covers  the  modelling  of  a  simple  autopilot,  which  is  important 
for  simulating  lateral  plane  motions  in  waves.  Section  7  describes  the  numerical 
implementation  of  the  present  work,  which  has  been  validated  using  experimental 
data  for  a  steered  warship  model  in  Section  8.  Recommendations  for  future  work 
are  given  in  Section  9,  followed  by  final  conclusions  in  Section  10. 

2  Translating  Earth  Axis  System  for 
Computing  Forces  on  a  Hull  and 
Appendages 


As  was  done  in  previous  work  on  the  ShipMo3D  library  [1,  2,  3],  the  present  report 
uses  a  translating  earth  axis  system  that  moves  with  the  mean  steady  ship  speed. 
Figure  1  shows  the  translating  coordinate  system.  The  relative  sea  direction  (3  uses 
the  convention  shown  in  Figure  2,  with  f3  =  180°  representing  head  seas. 

Note  that  the  translating  earth  axis  system  in  the  present  work  differs  from  the 
stability-based  axis  system  used  by  Schmitke  [4]  for  deriving  forces  on  a  ship  and 
its  appendages.  The  stability-based  axis  system  rotates  with  the  instantaneous  yaw 
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motion  of  the  ship;  thus,  yaw  displacement  and  forces  dependent  on  yaw  displace¬ 
ment  are  always  zero  when  using  a  stability-based  axis  system.  In  contrast,  forces 
in  a  translating  earth  axis  system  include  terms  proportional  to  yaw  displacement, 
such  as  lift  stiffness  terms  due  to  the  angle  of  attack  caused  by  yaw  displacement. 

3  Forces  on  Appendages 


Forces  on  appendages  significantly  influence  lateral  plane  motions  of  ships.  Much 
of  the  present  work  is  based  on  Schmitke  [4].  In  this  section,  radiation  and  dif¬ 
fraction  forces  due  to  appendages  are  assumed  to  be  negligible.  Computations  by 
McTaggart  and  Stredulinsky  [6]  for  HMCS  NIPIGON  confirmed  that  radiation  and 
diffraction  forces  on  appendages  will  have  negligible  influence  on  predicted  mo¬ 
tions  for  naval  destroyers.  The  evaluation  of  radiation  and  diffraction  forces  on 
appendages  is  discussed  in  Section  5  of  this  report. 


3.1  Foils  with  Small  Chord  Lengths  -  Rudders, 
Propeller  Shaft  Brackets,  and  Stabilizer  Fins 

Appendages  such  as  rudders,  propeller  shaft  brackets,  and  stabilizer  fins  are  as¬ 
sumed  to  have  small  chord  lengths.  Figures  3  and  4  give  characteristic  dimensions 
for  foils  with  small  chord  lengths. 

For  predicting  ship  motions  in  waves,  a  foil  chord  is  considered  to  be  small  if  it 
is  small  relative  to  both  the  length  of  the  ship  and  to  incident  wavelengths.  The 
assumption  of  a  small  chord  length  leads  to  the  following: 

•  lift  forces  acting  on  the  foil  motion  act  at  distance  c/4  forward  of  the  foil  cen¬ 
troid,  where  c  is  the  mean  chord  length, 

•  lift  forces  due  to  foil  motion  are  based  on  the  motion  at  distance  c/4  aft  of  the 
foil  centroid, 

•  forces  due  to  foil  added  mass  are  based  on  the  acceleration  of  the  foil  centroid, 

•  forces  due  to  foil  added  mass  act  at  the  foil  centroid. 

The  added  mass  of  a  foil  moving  perpendicular  to  its  plane  is  follows: 

AfpM  =  p\  sc2  (1) 


2 
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where  p  is  water  density.  The  lift  force  acting  on  a  foil  with  flow  velocity  V  at 
incident  angle  a  is  evaluated  by: 


rpfoil 
”  lift 


\pSV2a 


dClift 

da 


(2) 


where  S  is  the  foil  area  and  dCh^ /da  is  the  lift  curve  slope.  Equation  (2)  excludes 
a  circulation  delay  function  from  Schmitke  [4].  The  circulation  delay  function  is 
not  included  in  the  ShipMo3D  library  because  its  value  is  typically  close  to  one 
and  because  it  is  difficult  to  include  it  consistently  in  frequency  domain  and  time 
domain  formulations.  For  lift  forces  acting  on  a  ship  foil,  the  following  relationship 
is  used  for  evaluating  the  incident  flow  speed  and  angle: 


V2  a 


U 


A  Ufoil\ 

~U~  ) 


Vp 


(3) 


where  U  is  ship  forward  speed,  AU^ml  is  the  local  velocity  correction  at  the  foil, 
and  Vp  is  flow  velocity  perpendicular  to  the  foil. 


The  lift  curve  slope  in  uniform  flow  can  be  evaluated  using  the  following  equation 
from  Whicker  and  Fehlner  [7]: 


dChft 

da 


0.9  x  4w  s/c 
1.8  +  yj  (2s/c)2  +  4.0 


(4) 


Figure  5  shows  the  variation  of  lift  coefficient  slope  with  aspect  ratio  using  the 
above  equation,  and  also  the  lift  coefficient  slope  based  on  low  aspect  ratio  theory 
as  follows: 


dCHft 

da 


7T  s/c 


(5) 


Using  the  terms  given  above,  the  lift  damping  due  to  a  foil  moving  with  velocity 
perpendicular  to  its  plane  can  be  expressed  as: 


Tjfoil 

p 


A  Ufoil\ 

~U~  ) 


dClift 

da 


(6) 


The  local  incident  velocity  ratio  (l  +  AU^ml/U)  will  often  be  less  than  1 .0  due  to 
the  influence  of  the  hull  boundary  layer.  Propeller-induced  flow  and  the  influence 
of  the  hull  on  the  incident  potential  flow  will  also  influence  the  velocity  increment 
AUf°d.  SHIPM07  [8]  includes  the  effect  of  the  hull  boundary  layer  when  evalu¬ 
ating  local  incident  flow  velocity.  The  turbulent  boundary  layer  can  be  modelled 
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using  a  one-seventh  power  approximation  [9],  with  the  boundary  layer  thickness 
given  by: 


6b 


0.377  (: xFP  -  xfoil ) 
(Ref™1)0'2 


(7) 


The  foil  Reynolds  number  is  dependent  on  the  distance  of  the  foil  from  the  hull 
forward  perpendicular  as  follows: 


,  U  (xFP  -  xfoil) 

Refml  =  — ^-771 - L  (8) 

v 

where  xFp  is  the  longitudinal  location  of  the  forward  perpendicular,  xjmi  is  the 
longitudinal  location  of  the  foil,  and  v  is  the  kinematic  viscosity  of  water.  The 
speed  correction  factor  to  account  for  the  influence  of  the  hull  boundary  layer  is 
given  by: 


(1  +  AUfoil/U) 


1  + 


2  sc, 


root 


8  s 


Ctip) 


4  Croot 


(9) 


Figure  6  shows  the  influence  of  the  hull  boundary  layer  on  lift  of  rectangular  foils 
near  the  stem  of  a  frigate.  As  ship  speed  increases,  the  boundary  layer  thickness 
and  corresponding  speed  reduction  decrease.  As  foil  span  increases,  the  relative 
importance  of  the  boundary  layer  decreases. 
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Ship  Speed  (knots) 


Figure  6:  Influence  of  Hull  Boundary  Layer  on  Lift  Forces  from  Rectangular  Foils 
Near  Stern  of  Frigate  -  100  m  aft  of  Fore  Perpendicular 


For  roll  motions,  ship  damping  forces  due  to  potential  flow  are  often  small  and 
viscous  forces  can  be  significant.  Viscous  forces  on  foils  are  evaluated  by  treating 
them  as  flat  plates  oriented  at  angles  of  high  incidence  to  flow,  with  drag  forces  as 
follows: 

F-f°u~msc  =  l-pV*scCd  (10) 

where  Cd  is  the  drag  coefficient,  which  is  assumed  to  be  1.17  for  a  flat  plate. 

Using  the  forces  described  above,  hydrodynamic  coefficients  for  foils  appended  to 
ships  can  be  evaluated.  Each  foil  has  a  dihedral  angle  Y  relative  to  the  y  axis  as 
shown  in  Figure  4.  Each  foil  has  an  effective  roll  radius  for  lift  forces  based  on  the 
location  of  its  centroid  given  as  follows: 

ri  =  |7  cos  T  +  Tsinr  (11) 

For  roll-roll  added  mass  and  damping,  the  nominal  square  value  of  the  effective  roll 
radius  is  required,  and  is  evaluated  as  follows: 

rf  =  4  /  ( ycosY  +  zsinT)2  dS  (12) 

Js 

The  resulting  hydrodynamic  coefficients  in  translating  earth  axes  are: 

A[f  =  Afpoil  sin2  T  (13) 
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to  foil 
D22 

= 

Bfpoil  sin2  T 

(14) 

a  foil 
^23 

= 

—Afpml  cosT  sinT 

(15) 

jd  foil 
n23 

= 

-Bfpoil  cos  T  sinT 

(16) 

A  foil 
^24 

= 

-A**  sin  T  rj 

(17) 

jofoil 

n2A 

= 

-Bfpoil  sin  r  r, 

(18) 

A  foil 
^25 

= 

-TAfoil 

x/i23 

(19) 

to  foil 
■d25 

= 

-  (*  -  5/4)  S4“‘ 

(20) 

so  foil 
°25 

= 

f  B'f 

(21) 

A  foi1 
^26 

= 

t  Afoil 

(22) 

to  foil 
-°26 

= 

(X  -  c/4) 

(23) 

/-of  oil 
l'26 

= 

-U 

(24) 

A  foil 

^32 

= 

cosT  sinT 

(25) 

to  foil 
n32 

= 

-B^1  cosT  sinT 

(26) 

A  foil 
^33 

= 

Af*1  cos2r 

(27) 

TO  foil 
n33 

= 

Bspml  cos2r 

(28) 

A  foil 
^34 

= 

A^ml  COS  r  T'l 

(29) 

to  foil 
-°34 

= 

Bfpoil  cos  r  n 

(30) 

A  f°U 

^35 

= 

-T  Afoil 

x  ^1 33 

(31) 

jo  foil 
n35 

= 

-  (3T  -  c/4)  Bg" 

(32) 

so  foil 
°35 

= 

U  B{f 

(33) 

A  foil 

^36 

= 

t  Afoil 
x  4I32 

(34) 

jo  foil 
n36 

= 

(*  -  c/4)  B/f 

(35) 

so  foil 
°36 

= 

-V  B(f 

(36) 

A  foil 

^42 

= 

-Afpoil  sin  Yn 

(37) 

to  foil 
-°42 

= 

-Bfpoil  sin  T  r, 

(38) 

A  foil 

^43 

= 

A^ml  COS  r  T'l 

(39) 

to  foil 
nA3 

= 

Bfpoil  cos  rr, 

(40) 

A  foil 
41 44 

= 

r2 

(41) 

to  foil 
-D44 

= 

Bfp°il  rf  + 

(42) 

A  foil 

^45 

= 

X  /I43 

(43) 

to  foil 
-°45 

= 

-  -  c/4)  54“' 

(44) 

so  foil 
°45 

= 

U  B{f 

(45) 
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A  foil 
^46 

Tf.  Afoil 

—  X  -TI42 

(46) 

id  foil 
-°46 

=  (X  -  c/4)  B{f 

(47) 

nfoil 

*~'46 

=  -U  B{f 

(48) 

A  foil 

-  t  Afoil 

(49) 

td  foil 
-°62 

=  (x  +  c/2)  Bf22l 

(50) 

a  foil 
^63 

T  Afoil 
~  X  ^23 

(51) 

TDfoil 

-°63 

=  (x  +  c/2)  Bf2°d 

(52) 

a  foil 
^64 

-  -f  Afoil 

—  A  ^*-24 

(53) 

jo  foil 
-°64 

=  (x  +  c/2)  Bf2f 

(54) 

a  foil 
^65 

-gi  Afoil 
~  x  ^23 

(55) 

to  foil 
n65 

=  ~  (x  +  c/2)  (x+c/2)  B{ f 

(56) 

r>f°il 

°65 

=  UB/f 

(57) 

a  foil 
^66 

—  7vt2  Afoil 

—  x  +i22 

(58) 

TD  foil 

nm 

=  (x  -  c/4)  ( x  +  c/4)  B+ 

(59) 

si  foil 
°66 

=  -UB>? 

(60) 

The  roll  damping  term  B{^1  includes  a  contribution  B{™1  vlsc  from  viscous  damp¬ 
ing  which  can  be  evaluated  based  on  Equation  (10).  The  viscous  roll  moment  in¬ 
duced  by  roll  velocity  r) 4  will  be: 

F/“'-~c  =  \p  ills  Ctrl  (61) 

where  is  the  effective  cube  of  the  roll  moment  arm  for  viscous  forces  given  by: 

rl  —  —  f  (y  cosT  +  z  sinT)2y/?/2  +  z2  dS  (62) 

<5  Js 

For  computations  in  the  time  domain,  the  viscous  roll  damping  coefficient  is: 

B{f~msc  =  1  p  rh  S  Cd  rl  (63) 

For  computations  in  the  frequency  domain,  the  equivalent  linear  damping  coeffi¬ 
cient  based  on  energy  dissipation  is: 

Bfoii-visc  =  4  p  UJe  fj4  r3  g  Qd  (64) 

D7T 

where  f/4  is  the  roll  amplitude. 

When  foils  are  in  the  presence  of  waves,  excitation  forces  will  occur  due  to  added 
mass  and  lift  in  the  presence  of  incident  flow.  The  excitation  forces  can  be  computed 
using  the  added  mass  and  damping  terms  previously  computed: 

Ffoil  =  ay  +  az  A^f  +  uy  Bf™1  +  uz  Bjf  (65) 
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where  ay  and  az  are  the  horizontal  and  vertical  wave-induced  flow  accelerations, 
and  uy  and  uz  are  the  horizontal  and  vertical  wave-induced  flow  velocities. 


Some  foils  such  as  rudders  and  stabilizing  fins  will  have  deflections  which  will 
induce  forces  on  the  ship.  Schmitke  [10]  presents  forces  arising  from  rudders  and 
stabilizing  fins.  In  the  current  work,  forces  proportional  to  deflection  acceleration 
and  velocity  are  assumed  to  be  negligible.  The  convention  used  for  foil  deflection 
is  positive  for  counter-clockwise  when  viewed  from  inside  the  ship.  For  example, 
a  conventional  vertical  rudder  will  have  a  positive  deflection  angle  if  it  is  counter¬ 
clockwise  when  viewed  from  above.  The  lift  force  acting  perpendicular  to  a  rudder 
due  to  its  deflection  §rudder  will  be: 


7! rudder 
lift 


1 

_  _  q  ta2  rrudder 

2f  da 


The  stiffness  coefficient  for  the  rudder  deflection  can  be  expressed  as: 


(jrudder 


=  2»SU2 


A  Ufoil\2  dClift 
1  +  U  )  da 


Stiffness  coefficients  relating  ship  forces  to  rudder  deflections  are: 


s~i rudder 
°2  5 
irudder 
°3  5 
rudder 
°4  5 

s~irudder 
°5  5 

C  rudder 
65 


Cr5udder  sin  T 

-crsudder  cos  r 

s~i rudder 
—  L/<5  rl 

-Cl  (x  +  c/2) 
Cr2fder  (x  +  c/2) 


(66) 


(67) 


(68) 

(69) 

(70) 

(71) 

(72) 


3.2  Bilge  Keels  and  Other  Foils  with  Long  Chord 
Lengths  and  Low  Aspect  Ratios 

So  far  discussion  has  been  limited  to  foils  with  small  chord  lengths.  In  contrast, 
bilge  keels  have  long  chord  lengths,  which  are  non-negligible  relative  to  the  length 
of  the  ship  and  incident  wavelengths.  Bilge  keels  are  assumed  to  have  low  aspect 
ratios  with  lift  curve  slope  given  by  Equation  (5).  Figure  5  suggests  that  the  assump¬ 
tion  of  low  aspect  ratio  is  valid  for  aspect  ratios  less  than  0.3,  which  is  consistent 
with  typical  bilge  keels.  The  present  discussion  is  also  applicable  to  skegs  and  other 
appendages  with  long  chord  lengths  and  low  aspect  ratios. 

The  added  mass  of  a  bilge  keel  moving  perpendicular  to  its  plane  is  determined 
from  the  following  equation  based  on  potential  flow  theory  for  a  flat  plate  [11]: 

_  f*  f  OT*G 

Af  =  p-  s2(x)dx  (73) 

^  dxaft 
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where  xaft  is  the  aft  extent  of  the  bilge  keel,  x jore  is  the  forward  extent  of  the  bilge 
keel,  and  s(x)  is  the  local  span. 


Newman  [9]  and  Crane,  Eda,  and  Landsburg  [12]  discuss  lift  forces  on  low  aspect 
ratio  foils.  Equation  (5)  indicates  that  the  lift  coefficient  slope  of  a  low  aspect  foil 
is  proportional  to  the  aspect  ratio.  Consequently,  the  lift  force  acting  on  a  bilge  keel 
will  be  independent  of  chord  length,  and  can  be  expressed  as: 

Fujt  =  \  P  sL,  V2  <*  tt  (74) 

If  the  bilge  keel  has  constant  span  (i.e.,  independent  of  x),  then  the  lift  force  will  be 
concentrated  at  the  leading  edge.  Note  that  the  lift  magnitude  and  effective  location 
will  be  unaffected  if  local  span  decreases  downstream  of  the  location  of  maximum 
span.  To  simplify  the  analysis  of  lift  forces  on  bilge  keels,  the  lift  force  is  evaluated 
at  the  foremost  location  at  which  the  maximum  span  occurs,  denoted  xuft. 

The  lift  damping  due  to  a  bilge  keel  moving  perpendicular  to  its  plane  can  be  ex¬ 
pressed  as: 

B?  =  \  P  V  (1  +  A£/“/£7)  7,  (75) 

where  AUbk  is  the  local  incident  velocity  correction  at  the  bilge  keel. 

Viscous  eddy  damping  from  bilge  keels  can  be  a  significant  fraction  of  total  roll 
damping,  particularly  at  low  ship  speeds.  The  viscous  drag  force  on  a  bilge  keel 
can  be  expressed  as  follows: 

Fbk~visc  =  \pCd  f  xforeVp\x)s\x)dx  (76) 

^  J  Xaft 


Forces  acting  on  bilge  keels  are  significantly  influenced  by  two  flow  phenomena 
that  are  typically  neglected  for  shorter  foils  of  larger  aspect  ratios.  Drag  coefficients 
for  oscillating  bilge  keels  are  often  much  greater  than  the  value  of  1.17  for  a  flat 
plate  in  uniform  flow.  The  higher  drag  coefficients  are  due  to  the  oscillating  nature 
of  the  flow,  and  increase  with  decreasing  ratio  of  foil  span  to  motion  amplitude. 
The  second  phenomenon  that  should  be  considered  is  the  influence  of  the  hull  on 
the  incident  flow  normal  to  the  bilge  keel  and  resulting  drag.  This  effect  can  be 
modelled  using  a  local  roll  velocity  ratio  v,  which  is  the  ratio  of  the  local  normal 
velocity  induced  by  roll  to  the  velocity  in  the  absence  of  the  hull. 

Using  the  forces  described  above,  hydrodynamic  coefficients  for  bilge  keels  can 
be  evaluated.  The  local  span  of  a  bilge  keel  is  assumed  to  be  small  relative  to  the 
distance  from  the  ship  center  of  gravity;  thus,  the  effective  y  and  z  locations  are 
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taken  at  mid-span  of  the  bilge  keel.  The  resulting  hydrodynamic  coefficients  in 
translating  earth  axes  are: 


7T 

P  2  I 

PX  fore 

j  s2(x)  sin2  T(x)  dx 

(77) 

xaft 

Bbk  sin2  T(xuft) 

(78) 

7 r 

- P  o 

XX  fore 

/  s2(a;)  sin  T(x)  cosT  (a;)  dx 

(79) 

Z 

dxaft 

-Bbk 

sin  T(xuft)  cos  T(x uft) 

(80) 

7 r 

-P  p 

XX  fore 

/  s2(a;)  sinT(a;)  ri{x)  dx 

(81) 

Z 

dxaft 

-Bbk 

sin  T(xuft)  r i(xuft) 

(82) 

P  —  s2(x)  x  sin  T(x)  cosT  (a;)  dx 

^  Jxaft 

Bbk  x  sin  F(xuft)  cos  T(xuft) 


s2(x )  sin2  T(x)  x  dx 


Xu  ft  Bbk 
-U  Bbk 

7j-  rxfore 

-P  o  / 


s2(x )  sinr(x)  cosT  (a;)  dx 


-Bbpk  smr(xuft)  cos T(xuft) 


p  —  /  s2(x)  coszT(x)  dx 

^  Jxaft 

Bbpk  cos2  T(xuft) 

fXfore 

„  "  /  „2/  \ 


s2(x)  cosT  r;(x)  da; 


Bp  cosT (x Hft)  ri(xHft) 


s2(x )  cosT  r;(x)  dx 


Bbk  cosT (x uft)  n(xuft) 

f*  f  ovg. 

p  -  s2(x)  ri(x)2  dx 
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Tjbk 

baa 

=  Bf  r?(xK„)  + 

(100) 

Abk 

^45 

7|-  Before 

=  —p  —  /  s2(x)  x  cosT  riix)  dx 

(101) 

^  Jxaft 

Tjbk 

-^45 

=  Bbpk  xuft  cosT(xiift)  ri(xuft) 

(102) 

/" <bk 
*~'45 

=  UBbk 

(103) 

Abk 

^46 

7j-  f‘xfore 

=  —  p  —  /  s2(x)  sin  T(x)  riix)  x  dx 

(104) 

^  Jxaft 

Tjbk 

-°46 

=  xliftBbk 

(105) 

/~<bk 

°46 

Abk 

^62 

=  - UBbk 

7j-  rxfore 

=  p—  /  s2(x)  sin2  T(x)  x  dx 

(106) 

(107) 

^  Ao/t 

rybk 

-^62 

= 

(108) 

Abk 

^63 

yj-  f‘x  fore 

=  —  p  —  /  s2(x)  sin  T(x)  cos  T (a;)  a;  dx 

(109) 

^  J  Xaft 

r>bk 

-^63 

=  *uftBbk 

(HO) 

Abk 

^64 

fx  fore 

=  —  p  —  /  s2(x)  sin r(x)  riix)  x  dx 

(111) 

^  J  xaft 

Tjbk 

-^64 

=  xliftBbkA 

(112) 

_  f*  X  j-  OVG 

^65  =  p—  /  s2(x)  sinr(x)  cos  F(x)  x2  dx  (113) 

2  ^ 


xaft 


r>bk 

^65 

■r2  Rbk 

~  Xh ft  n23 

(114) 

/~<bk 

065 

Abk 

^66 

=  UBbk 

7j-  rxf  ore 

=  p—  /  s2(x)  sin2  T(x)  a:2  dx 

(115) 

(116) 

^  Jxaft 

rybk 

-^66 

2  n&fc 

—  xZi/t  -°22 

(117) 

/^ibk 

°66 

=  U  Bq2 

(118) 

The  effective  square  of  the  local  added  mass  moment  arm  is  given  by: 

rf(x)  =  [ V(x )  cosT(a;)  +  z(x)  sinT(a;)]2  (119) 

The  roll  damping  term  B\\  includes  a  contribution  BbA4~msc  from  viscous  damping 
which  can  be  evaluated  based  on  Equation  (76).  The  viscous  roll  moment  induced 
by  roll  velocity  /)4  will  be: 

1  P^fore  ^ 

pbk-viBc  =  pf)lv2  /  s(x)  Cd(x)  r3(x)  dx  (120) 

^  JXaft 
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where  r%(x)  is  the  effective  cube  of  the  local  roll  moment  arm  given  by: 


r%(x)  =  \y(x)  cosr(a;)  +  z(x)  sinr(a;)]2  \Jy2(x )  +  z2{x)  (121) 

For  computations  in  the  time  domain,  the  viscous  roll  damping  coefficient  is: 

1  fXfore  ^ 

Bbtvisc  =  -PV 4  /  s(x)  Cd(x)  r*(x)  dx  (122) 

For  computations  in  the  frequency  domain,  the  equivalent  linear  damping  coeffi¬ 
cient  based  on  energy  dissipation  is: 

A  r^fore  , _ , 

Bbu~msc  =  —puef)4  /  s(x)  Cd(x)  r%(x)  dx  (123) 

3?r  Jxaft 


Excitation  forces  acting  on  bilge  keels  arise  from  both  added  mass  and  lift  compo¬ 
nents,  and  are  evaluated  as  follows: 

Ff  =  ay  Af2  +  az  A%  +  uy  B%  +  uz  Bf3  for  j  =  2  -  6  (124) 

The  flow  velocity  components  CX/y  j  CX/ g  j  XLy  j  and  uz  are  evaluated  at  the  bilge  keel 
centroid  of  area.  For  appendages  which  are  long  relative  to  incident  wavelengths,  it 
is  recommended  that  that  excitation  forces  be  evaluated  using  several  longitudinal 
segments. 

3.3  Drag  Coefficients  and  Roll  Velocity  Ratios  for 
Bilge  Keels 

One  of  the  most  challenging  problems  in  seakeeping  computations  is  the  selection 
of  appropriate  drag  coefficients  for  bilge  keels.  Note  that  the  present  discussion  is 
also  applicable  to  other  long  foils  of  low  aspect  ratio,  such  as  skegs.  Roll  motion 
predictions  can  be  very  sensitive  to  input  drag  coefficients  for  bilge  keels.  Further¬ 
more,  different  prediction  methods  often  give  widely  varying  results  for  bilge  keel 
drag  coefficients.  One  common  feature  among  results  from  different  prediction 
methods  is  that  the  drag  coefficient  increases  with  decreasing  roll  amplitude.  At 
small  roll  amplitudes,  experimental  results  [5]  indicate  that  bilge  keel  drag  coeffi¬ 
cients  will  often  exceed  10,  a  surprising  finding  considering  that  the  drag  coefficient 
of  a  flat  plate  normal  to  uniform  flow  is  approximately  1.17.  The  high  observed 
drag  coefficients  for  bilge  keels  are  due  to  the  oscillatory  nature  of  the  bilge  keel 
motions,  with  drag  coefficients  generally  increasing  as  the  ratio  of  the  bilge  keel 
span  to  motion  amplitude  decreases. 
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Ikeda’s  method  for  predicting  bilge  keel  drag  coefficients  [5,  13]  is  widely  used. 
The  drag  coefficient  is  given  by: 


Cd 


22.5 - - -  +  2.4 

7T  rBK  f)4 


(125) 


where  rBK  is  the  roll  radius  from  the  ship  centre  of  gravity  to  the  bilge  keel.  Ikeda’s 
method  uses  the  following  equation  for  the  the  local  roll  velocity  ratio: 


v 


1.0  +  0.3  exp 


-160  1.0 


Ax  \ 
B  x  Tx) 


(126) 


where  Ax  is  hull  cross-sectional  area,  Bx  is  hull  sectional  beam,  and  Tx  is  hull 
sectional  draft.  Figure  7  shows  the  local  roll  velocity  ratio  as  a  function  of  hull 
sectional  area  coefficient. 


SHIPM07  [8]  and  it’s  predecessors  use  a  method  developed  by  Kato  [14]  and  de¬ 
scribed  by  Schmitke  [4]  for  computing  bilge  keel  drag  coefficients.  This  method 
appears  to  be  overly  complex  given  the  relatively  limited  number  of  parameters 
upon  which  it  is  based.  Using  calculations  for  the  destroyer  HMCS  NIPIGON  from 
McTaggart  and  Stredulinsky  [6],  the  following  simplified  equations  can  be  devel¬ 
oped  from  Kato’s  method: 

Cd  =  5  (127) 

V  =  J  1.0  +  3.5  exp  ^  bilge^j  (128) 


where  Rbuge  is  the  hull  bilge  radius.  In  contrast  to  Ikeda’s  method,  Kato’s  equation 
for  drag  coefficient  has  a  roll  velocity  amplitude  term  //4  instead  of  a  roll  amplitude 
term  //4.  Furthermore,  Kato’s  method  includes  gravitational  acceleration  g.  Due 
to  the  submergence  of  typical  bilge  keels  and  resulting  lack  of  radiated  waves,  it 
is  questionable  whether  gravitational  acceleration  actually  influences  their  viscous 
roll  damping.  However,  good  roll  predictions  for  HMCS  NIPIGON  [6,  15]  ob¬ 
tained  using  Kato’s  formulation  suggest  that  the  approach  has  merit.  Given  that 
Kato’s  method  gives  good  results  for  HMCS  NIPIGON,  with  bilge  keel  parameters 
given  in  Table  1,  the  following  new  equation  is  proposed  for  the  drag  coefficient 
dependent  on  roll  amplitude  rather  than  roll  velocity  amplitude: 

Cd  =  16  (129) 

It  should  be  emphasized  that  Equations  (127)  and  (129)  are  initial  attempts  at  sim¬ 
plified  approaches  based  on  Kato’s  method.  Ideally,  numerical  values  in  Equations 
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Table  1:  Bilge  Keel  Parameters  for  HMCS  Nipigon 


Length  c  26.5  m 

Span  s  0.61  m 

Average  bilge  keel  radius  fBi<  6.8  m 
Average  bilge  hull  radius  Ruige  3.2  m 
Ship  natural  roll  period  10.6  s 


(127)  and  (129)  should  be  optimized  based  on  analysis  of  large  volumes  of  experi¬ 
mental  data. 


To  compare  local  roll  velocity  ratios  from  Ikeda’s  and  Kato’s  methods,  it  is  useful 
to  consider  rectangular  hull  sections  with  rounded  corners.  The  cross-sectional  area 
is  given  by  the  following: 

A,  =  Bx  T„  -  (2  -  |)  R2Wge  (130) 

Figures  7  and  8  show  local  roll  velocity  ratios  as  functions  of  sectional  area  co¬ 
efficient  and  nondimensional  bilge  radius  for  rounded  rectangular  sections  with  a 
beam  Dx  of  8  m,  a  draft  Tx  of  4  m,  and  the  center  of  gravity  located  1  m  above 
the  waterline.  Bilge  keel  radii  tbk  are  evaluated  at  the  mid-span  location  for  bilge 
keels  with  spans  of  0.6  m.  For  sections  with  less  rounded  corners  (i.e.,  larger  area 
coefficients  and  smaller  non-dimensional  bilge  radii),  local  roll  velocity  ratios  from 
Kato’s  method  are  much  greater  than  those  from  Ikeda’s  method.  Fortunately,  the 
influence  of  differences  between  the  two  methods  is  mitigated  by  the  fact  that  bilge 
keel  roll  damping  is  typically  of  lesser  relative  importance  for  hull  sections  with 
small  bilge  radii. 


As  an  alternative  to  the  methods  of  Ikeda  and  Kato,  Lloyd  [16]  presents  plotted 
bilge  keel  drag  coefficients  and  associated  equations  based  on  experimental  data 
Unfortunately,  Lloyd’s  equations  appear  to  have  errors  while  the  plotted  results 
appear  to  be  correct.  The  following  equations  have  been  developed  to  represent 
Lloyd’s  plotted  results: 


Cd 

dCd 

bed 


1.17+11.23  acd 


V 4  tBk 


11.23  exp 


0.3  yfsc \ 

0.0125  +  0.025  774  tbk/Vs~c 


(131) 

(132) 

(133) 


For  Equations  (132)  and  (133),  the  roll  amplitude  r/4  is  given  in  radians. 
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Area  Coefficient  AX/(BX  Tx  ) 


Figure  7:  Local  Roll  Velocity  Ratio  Versus  Sectional  Area  Coefficient  for 
Rectangular  Hull  Section  with  Rounded  Corners 


Nondimensional  Hull  Bilge  Radius  Rbiige/rBK 


Figure  8:  Local  Roll  Velocity  Ratio  Versus  Nondimensional  Hull  Bilge  Radius  for 
Rectangular  Hull  Section  with  Rounded  Corners 
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Figure  9:  Bilge  Keel  Drag  Coefficients  based  on  Lloyd  [16] 


Figure  10  shows  predicted  bilge  keel  drag  coefficients  versus  roll  amplitude  for 
HMCS  NIPIGON.  There  are  significant  differences  among  the  3  different  methods, 
particulary  at  low  roll  amplitudes  for  which  Lloyd’s  method  gives  much  smaller 
drag  coefficients. 

4  Viscous  Hull  Roll  Damping 


When  a  ship  hull  rolls  in  waves,  significant  damping  forces  arise  from  eddy  mak¬ 
ing  and  skin  friction.  Schmitke  [4]  and  Graham  [17]  discuss  prediction  of  hull 
eddy  making  and  skin  friction  forces  within  strip  theory  computations.  This  section 
presents  methods  for  computing  viscous  hull  roll  damping  for  a  panelled  ship  hull. 

4.1  Formulation  for  Hull  Eddy-Making  Roll  Damping 

SHIPM07  [8]  and  its  predecessors  evaluate  hull  eddy-making  roll  damping  using 
Tanaka’s  method  [18].  Graham  [17]  gives  a  detailed  description  of  the  SHIPMO 
implementation.  Tanaka’s  method  is  based  on  empirical  expressions  based  on  sec¬ 
tional  geometric  parameters. 

A  new  approach  is  introduced  here  for  computing  eddy-making  roll  damping  using 
a  hull  geometry  described  by  surface  panels.  Like  Tanaka’s  method,  the  new  ap¬ 
proach  correctly  models  the  limiting  cases  of  a  circular  cylindrical  hull  with  zero 
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Figure  10:  Predicted  Bilge  Keel  Drag  Coefficient  Versus  Roll  Amplitude  for  HMCS 
NIPIGON 

eddy-making  damping  and  a  hull  acting  as  a  flat  plate.  Using  a  flat  plate  analogy, 
the  eddy-making  damping  moment  at  zero  ship  speed  is  given  by: 

FZ'u‘-*dd>  =  \p  |,j4|  CX  l  [  V?  +  *2  dS  (134) 
2  2  JshM 

where  C^dy  is  the  hull  drag  force  coefficient  and  n4  is  the  local  roll  normal  given 
by: 


n4  =  y  nz  —  z  ny  (135) 

As  a  first  approximation,  the  drag  coefficient  can  be  given  a  value  of  1.17 
based  on  a  flat  plate  moving  normal  to  flow.  The  fraction  of  1/2  preceeding  the 
integral  in  Equation  (134)  accounts  for  the  integral  including  both  the  starboard 
and  port  sides  of  the  hull. 

For  computations  in  the  time  domain,  the  hull  eddy-making  roll  damping  coeffi¬ 
cient  can  be  developed  from  Equation  (134).  Ikeda  et  al.  [19]  indicate  that  hull 
eddy-making  damping  decreases  significantly  with  forward  ship  speed,  and  have 
developed  the  following  relationship: 

Blf-eddv(U)  =  _ 1 _ 

Bhull-eddy( o)  1  +  625  iP/iLppUef 
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where  Lpp  is  ship  length  between  perpendiculars.  The  resulting  roll  damping  coef¬ 
ficient  including  speed  dependence  is: 


■pyhull — eddy  / t  t\  -i  n 

Bhu-My  =  Rt u-,mL[  I  P  fel  CS,  /  (137) 

For  computations  in  the  frequency  domain,  the  equivalent  linear  damping  coeffi¬ 
cient  based  on  energy  dissipation  is: 


B 


hull— eddy 
44 


t)  hull— eddy 

^44 _ 

r>  hull— eddy 
£>44 


SO  A 

(0)  3vr 


pcuf 


yhull 


V4  Ceddy 


z2  dS  (138) 


When  computing  ship  motions  in  the  time  domain,  the  wave  encounter  frequency 
will  often  be  varying  with  time.  It  is  suggested  that  the  ship  roll  natural  frequency 
be  used  for  the  speed  correction  term  in  Equation  (136)  because  roll  damping  is 
most  important  in  the  vicinity  of  roll  resonance. 


4.2  Drag  Coefficient  for  Hull  Eddy-Making  Roll 
Damping 

The  method  presented  above  requires  selection  of  suitable  value  for  the  hull  eddy 
damping  coefficient  C^ddy.  As  noted  above,  the  drag  coefficient  value  of  1.17  for  a 
flat  plate  is  suggested  as  an  initial  estimate  for  the  hull  eddy  damping  coefficient. 
Experimental  data  can  provide  further  insight  into  selection  of  suitable  values  for 
the  coefficient  C^ddy. 


4.3  Hull  Skin  Friction  Roll  Damping 

Additional  roll  damping  arises  from  flow  tangential  to  the  hull  surface.  The  roll 
moment  due  to  skin  friction  is  given  by: 

phuu-skin  =  I  p  \fjA\  fj4  C)uli  [  (\yn2\  +\zn3\  f  yV  +  z2  dS(  139) 

^  Jshull 

The  skin  friction  coefficient  is  calculated  using  the  Schoenherr  line  for  smooth  tur¬ 
bulent  flow: 

-2 

for  U  >0  (140) 

For  the  case  of  zero  forward  speed,  Schmitke  [4]  used  an  empirical  equation  from 
Kato  [20]  for  the  skin  friction  coefficient.  To  avoid  discontinuities  in  skin  friction 
coefficients  at  low  speeds,  a  new  approach  is  introduced  here  whereby  a  minimum 


C)ul1  =  0.0004  + 


3,46  log10  I  |  _  5.6 


V 
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ship  speed  threshold  is  introduced  for  Equation  (140).  This  approach  is  reasonable 
considering  the  slow  variation  of  skin  friction  coefficient  with  ship  speed.  For  a 
ship  length  of  100  m,  a  minimum  ship  speed  of  approximately  1  m/s  is  specified. 
Applying  Froude  scaling,  the  minimum  input  ship  speed  for  Equation  (140)  is: 

Urnvn(Chful1)  =  0.03  (141) 

The  resulting  roll  damping  coefficient  arising  from  skin  friction  is: 

Bkf~skin  =  \p  1 7)4 1  C)ul1  [  ( \yn2\  +  \z  n3|)2  ^Jy2  +  -2  dS  (142) 

^  J  shuH 

For  computations  in  the  frequency  domain,  the  equivalent  linear  damping  coeffi¬ 
cient  based  on  energy  dissipation  is: 

=  -i-  pcvefj4Ckul1  [  (; yn2  +  z  n3)2  y/y2  +  -2  dS  (143) 

d  71  Jshull 

5  Inclusion  of  Appendages  in  Hull 

Radiation  and  Diffraction  Calculations 


As  stated  in  Section  3,  ship  appendages  are  normally  assumed  to  have  negligible 
influence  on  ship  radiation  and  diffraction  forces.  This  assumption  greatly  sim¬ 
plifies  radiation  and  diffraction  calculations,  and  is  likely  valid  when  appendages 
are  deeply  submerged  and  very  small  relative  to  the  ship  hull.  In  efforts  to  ob¬ 
tain  improved  ship  motion  predictions,  a  method  has  been  developed  to  include 
appendages  in  hull  radiation  and  diffraction  calculations. 

Computations  for  HMCS  NIPIGON  [6]  showed  that  inclusion  of  appendages  in  hull 
radiation  and  diffraction  calculations  had  negligible  influence  on  predicted  motions. 
However,  appendage  radiation  and  diffraction  forces  could  be  significant  for  ships 
with  appendages  that  are  larger  or  closer  to  the  free  surface. 

5.1  Modelling  of  Appendages  Using  Dipole  Panels 

When  computing  hull  radiation  and  diffraction  forces  in  the  frequency  domain,  the 
ShipMo3D  library  represents  the  wetted  ship  hull  surface  using  a  series  of  panels 
[1].  Flow  boundary  conditions  on  the  hull  are  satisfied  by  solving  for  panel  source 
strengths.  Upon  initial  consideration,  appendages  could  be  included  in  the  calcu¬ 
lations  by  representing  each  appendage  using  source  panels.  Unfortunately,  this 
approach  is  not  practical  due  to  the  small  thickness  of  the  appendages  relative  to 
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practical  panel  sizes  (see  Sarpkaya  and  Isaacson  [11]).  Instead,  the  appendages  can 
be  assumed  to  be  very  thin,  and  represented  using  dipoles  representing  the  flow 
across  each  panel.  Sarpkaya  and  Isaacson  [11],  Chakrabarti  [21],  and  Meyerhoff 
[22]  provide  useful  references  for  computation  of  flows  using  dipole  panels. 

When  including  dipole  appendages  in  radiation  and  diffraction  computations,  the 
approach  for  an  unappended  ship  described  in  [1]  must  be  modified.  The  velocity 
potential  induced  by  a  dipole  panel  is  represented  by: 

0(f)  =  l-i  [  Gdipole(x,xd,ke)dS  (144) 

4vr  /o  , 

dipole 

where  x  is  the  field  point  location,  //  is  the  dipole  strength,  taken  as  being  constant 
over  the  panel,  Sdipoie  is  the  surface  of  the  dipole,  Gdipole  is  the  frequency  domain 
Green  function  for  a  dipole,  xd  is  the  location  on  the  dipole,  and  ke  is  the  wavenum¬ 
ber  based  on  encounter  frequency.  The  dipole  strength  //  is  equal  to  the  difference 
in  source  strength  between  the  front  and  back  of  a  dipole  panel.  The  Green  function 
for  a  dipole  is  related  to  the  Green  function  for  a  source  presented  in  Reference  1 
as  follows: 


Gdipole(x,xd,ke 


d G(x,  xd,  ke ) 

Qr^dipole 


(145) 


where  G(x ,  xd,  ke)  is  the  Green  function  for  a  source  and  ndipole  is  the  normal  vector 
from  the  dipole  panel.  The  following  equation  is  used  for  evaluating  the  normal 
derivative  of  the  source  Green  function: 


dG(x,xd,  fce) 

Qj^dipole 


d G(x,xd,  ke) 

o  dipole 

OTix 


+ 


d G(x,  x d,ke) 

r\  dipole 

UTly 


+ 


dG(x,  xd,  ke) 

r\  dipole 
OTlz 


(146) 


where  ndzpole,  ndipole,  and  nfpoh:  are  the  directional  components  of  the  normal  from 
the  dipole  panel.  There  are  two  possible  choices  for  the  normal  directional  from  a 
dipole  panel.  The  selected  normal  direction  will  determine  the  sign  of  the  dipole 
strength  fi  computed  for  the  panel. 


As  indicated  by  Meyeroff  [22],  the  dipole  strength  //  must  equal  zero  at  the  unbound 
edges  of  a  thin  appendage.  This  condition  typically  leads  to  high  gradients  in  dipole 
strength  near  the  edges;  thus,  it  is  often  appropriate  to  have  smaller  panel  sizes  near 
the  edges  of  an  appendage. 


As  discussed  by  Sarpkaya  and  Isaacson  [11],  the  only  net  pressure  on  a  dipole 
panel  arises  from  its  own  dipole  strength.  Other  dipoles  and  sources  induce  equal 
pressures  on  both  sides  of  the  dipole,  leading  to  a  net  pressure  of  zero  across  the 
dipole  panel. 
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5.2  Derivatives  of  the  Frequency  Domain  Green 
Function  for  a  Source  at  Zero  Forward  Speed 

As  indicated  above,  the  Green  function  for  a  dipole  is  based  on  the  normal  derivative 
of  the  Green  function  for  a  source.  Note  that  derivatives  are  required  for  the  Green 
function  of  a  dipole;  thus,  second  derivatives  are  required  for  the  Green  function  of 
a  source. 


The  ShipMo3D  library  uses  the  Telste  and  Noblesse  [23]  formulation  of  the  Green 
function  in  the  frequency  domain  at  zero  forward  speed: 

G(x,xs,ke)  =  4  +  4“  +  G0(x,xs,ke)  (147) 

it  111 

where  R  is  the  distance  from  the  field  point  x  to  the  source  at  xs,  R\  is  the  distance 
from  the  field  point  x  to  the  image  of  the  source  at  xs,  and  G0  is  the  frequency 
dependent  term  of  the  Green  function.  Garrison  [24]  discusses  the  Green  function 
at  the  zero  and  infinite  frequency  limits.  The  zero  frequency  (i.e.,  ke  =  0)  limit  is 
denoted  G0  and  given  by: 

G0(x,  xs,  ke)  =  4  +  4“  for  ke  =  0  (u;e  =  0)  (148) 

it  ill 


At  infinite  frequency  (i.e.,  ke  =  oo),  the  Green  function  is  denoted  G^  and  given 
by: 

Goo{x,X8,  fee)  =  4  -  7T  f°r  ke  =  OO  (ue  =  Oo)  (149) 

it  111 

Accordingly,  the  frequency  dependent  portion  Go  has  the  following  limits: 

G0(x,xs,  ke)  =  0  for  ke  =  0  (ue  =  0)  (150) 

~  2 

G0(x,xs,  ke)  =  -  —  for  ke  =  oo  (cce  =  oo)  (151) 

Ki 

Based  on  Telste  and  Noblesse,  the  frequency  dependent  term  is  expressed  as  fol¬ 
lows: 


Gq(x,Xs,  ke) 


2  ke 


Ro{h,  v) 


i  tv  J0(h)  exp(n) 


(152) 


where  Ro(h,  v)  is  a  function  to  be  evaluated  and  J0(/i)  is  the  Bessel  function  of  the 
first  kind  of  order  zero.  The  function  arguments  h  and  v  are  as  follows: 


Rxy 

=  \J(x  -  xs)2  +  {x  -  xs )2 

(153) 

h 

ke  Rxy 

(154) 

V 

=  ke  (z  +  zs) 

(155) 
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Thus,  Rxy  is  the  horizontal  distance  from  a  field  point  to  a  source,  h  is  the  non- 
dimensional  horizontal  distance,  and  v  is  the  non-dimensional  elevation  of  the  field 
point  relative  to  the  image  of  the  source. 


The  derivatives  of  the  Green  function  are  required  for  computing  flow  velocities 
induced  by  a  source.  The  zero  frequency  portion  of  the  Green  function  has  the 
following  derivatives, 


dG0  _  ~(x~xs)  (x-xs) 
dx  ~  R 3  R:\ 

dGo  =  -jy-ys)  _  {y  -  Vs) 
dy  R3  R:\ 

dG0  _  ~{z  -  zs)  _  (z  +  zs) 

dz  ~  R3  R\ 


(156) 

(157) 

(158) 


Teltse  and  Noblesse  give  the  following  equations  for  the  derivatives  of  the  fre¬ 
quency  dependent  portion  of  the  Green  function 


dGo 

9  Rxy 

:  -2/,- 

l  Ri(h,  v)  —  mJi{h)  exp(n) 

(159) 

dGo 

dRxy  dG 

(160) 

dx 

dx 

9Rxy 

dGo 

dRxy  dG 

(161) 

dy 

dy 

dRxy 

dGo 

dz 

=  2  k2e 

+  R0(h,v)  ITT  J0(h)  exp(n) 

ke  R\ 

(162) 

where  Ri(h,  v )  is  a  function  to  be  evaluated,  and  J\  (h)  is  the  Bessel  function  of  the 
first  kind  of  order  one.  The  derivatives  of  Rxy  with  respect  to  x  and  y  are: 


d  Rxy  (*T  Xs') 

dx  Rxy 

ORXy  {y  ys) 

dy  Rxy 


(163) 

(164) 


When  evaluating  the  flow  velocities  induced  by  a  dipole,  the  second  derivatives  of 
the  Green  function  are  required.  Relevant  terms  are: 


d2G0 

dRly 


- 2k 3 


dR\(h,v )  .  dJAh)  .  . 

~Sh~  - 1 1  — exp(,,) 


(165) 
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The  above  equations  utilize  the  following  relation  given  by  Telste  and  Noblesse: 

dR0(h,  v 


dh 


=  -  Ri(h,v) 


(176) 


The  derivatives  dR0{h,v)/dv  and  dRi(h,v)/dh  are  evaluated  using  finite  differ¬ 
ence  methods,  which  allow  utilization  of  the  fast  algorithms  of  Telste  and  Noblesse. 
The  Bessel  function  derivative  dJi(h)/dh  is  given  by: 


d.h(h) 

dh 


J0(h)  —  '—1^1  for  >  o 
h 

0.5  for  h  =  0 


(177) 

(178) 
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5.3  Algebraic  Solution  of  Hull  Source  and 
Appendage  Dipole  Strengths 

The  solution  of  hull  source  and  appendage  dipole  strengths  is  very  similar  to  the 
approach  used  for  a  hull  alone  represented  by  source  panels  in  Reference  1.  Lateral 
symmetry  is  utilized  to  reduce  the  required  computational  time.  Special  care  must 
be  taken  if  a  ship  includes  vertical  appendages  on  the  centreline,  such  as  a  skeg 
or  single  rudder.  For  longitudinal  modes,  the  dipole  strengths  on  centreline  ap¬ 
pendages  are  zero;  thus,  centreline  dipole  strength  terms  should  be  excluded  when 
solving  for  longitudinal  modes. 


5.4  Verification  of  Radiation  Computations  for 
Appended  Hull 

To  verify  radiation  computations  for  a  hull  with  appendages,  sway  added  mass  has 
been  evaluated  for  a  thin  box  with  length  of  20  m,  depth  of  10  m,  and  width  of 
0.4  m.  A  first  set  of  computations  models  the  box  as  a  hull  with  no  appendages.  A 
second  set  of  computations  models  the  upper  half  as  a  box  and  the  lower  half  as  a 
thin  plate.  Figure  11  shows  the  two  different  configurations,  denoted  as  “Box”  and 
“Box-plate”.  For  both  configurations,  there  are  20  panels  along  the  length  of  the 
body  and  20  panels  along  the  total  depth  of  the  body. 


Box 


Box-plate 


Figure  11:  Cross-Sections  of  Thin  Box  for  Verification  of  Radiation  Computations 
with  Appended  Hull 

The  computed  added  mass  and  damping  coefficients  in  Figures  12  and  13  show 
reasonable  agreement  considering  that  the  modelled  geometries  are  somewhat  dif- 


26 


DRDC  Atlantic  TM  2004-227 


ferent.  This  agreement  suggests  that  the  radiation  computations  using  dipole  ap¬ 
pendage  panels  have  been  implemented  correctly. 


6  Control  of  Rudder  Deflections  with  a 
Simple  Autopilot 


The  lateral  motions  of  a  ship  in  waves  will  be  influenced  by  rudder  deflections, 
which  are  often  driven  by  an  autopilot.  This  section  presents  a  relatively  simple 
autopilot,  and  provides  a  starting  point  for  simulation  of  rudder  motions.  Future 
studies  will  investigate  more  complex  helm  control  systems. 


In  response  to  ship  motions,  a  rudder  will  deflect  according  to  prescribed  motion 
gains  and  the  response  characteristics  of  the  rudder  system.  Equations  of  rudder 
motions  are  presented  in  various  references,  including  References  10  and  12.  The 
commanded  rudder  deflection  angle  determined  by  the  ship  motions  and  autopilot 
gains  is  as  follows: 

VS**  =  E  (*&  %  +  kh  * l i  +  Vi)  (179) 

3= 1 

where  k'^  is  the  acceleration  gain  for  motion  mode  j,  k'h]  is  the  velocity  gain,  and 
k$j  is  the  displacement  gain.  The  rudder  response  characteristics  are  modelled  as 
follows: 


£ rudder 


,  r>  /■  rudder 

+  ^  Qs  (^8 


£ rudder 


i  ,  ,2  rr udder 

+  u5  d 


uig  5rcudder  (180) 


where  (s  is  nondimensional  damping  response  constant,  and  u>sj  is  the  rudder  fre¬ 
quency  response  constant.  The  rudder  motions  in  the  time  domain  can  be  easily 
evaluated  as  follows: 


er  udder  ,  ,2  /  rr udder  rr  udder  \  q  /■  ,  ,  rr  udder  /ion 

o  =  UJS  {oc  -  d  )  -  2  QujS0  (181) 


Motions  in  the  frequency  domain  are  evaluated  using  the  following  coupled  equa¬ 
tion  of  ship  and  rudder  motions: 
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(182) 


The  matrices  and  vectors  in  the  above  equation  are  partitioned  into  ship  motion 
and  rudder  deflection  portions,  and  the  effective  mass  matrix  [Me]  includes  con¬ 
tributions  from  both  ship  mass  and  added  mass.  Terms  for  the  influence  of  rudder 
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Sway  Damping  B22/(pLT2Lue)  <S  Sway  Added  Mass  A22/{pLT2) 


‘  0  1  2  3  4 

Encounter  frequency  toe 


12:  Sway  Added  Mass  for  Thin  Box  and  Thin  Box  with  Appended  Plate 


‘  0  1  2  3  4 

Encounter  frequency  toe 


Figure  13:  Sway  Damping  for  Thin  Box  and  Thin  Box  with  Appended  Plate 
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deflections  on  ship  motions  were  previously  given  in  Section  3.  Other  terms  are  as 
follows: 


Mgj  —  ~^s  kgj 

(183) 

Mis  =  1 

(184) 

BSj  =  kvSj 

(185) 

B55  =  2  (s  (V5 

(186) 

Csj  =  — 

(187) 

Css  =  a; s 

(188) 

The  simulation  of  a  rudder  control  system  including  autopilot  requires  selection 
of  suitable  input  values.  For  a  US  Coast  Guard  cutter  representative  of  modern 
frigate  design.  Smith  [25]  indicates  a  maximum  rudder  deflection  of  35  degrees 
and  maximum  rudder  rate  of  3  degrees  per  second.  For  modelling  of  a  conventional 
downward  rudder  using  ShipMo3D,  yaw  gain  and  yaw  velocity  gain  will  typically 
have  values  greater  than  zero.  Values  reported  by  Smith  indicate  that  yaw  autopilot 
gains  fitted  to  observed  human  pilot  behaviour  will  have  a  wide  range  of  variation. 

7  Numerical  Implementation 


As  with  other  components  of  the  ShipMo3D  library,  the  present  work  has  been 
implemented  using  the  Python  programming  language  [26,  27].  Python  is  a  high 
level,  object-oriented  language  that  has  a  significant  user  base  among  the  scientific 
community.  The  ShipMo3D  library  uses  the  Numeric  Python  library  [28]  which 
provides  fast  computational  speed  for  numerical  operations.  The  scientific  Python 
module  SciPy  (www.scipy.org)  provides  additional  routines  such  as  fast  solution  of 
systems  of  linear  equations. 


7.1  Convergence  to  Correct  Roll  Amplitude  for 
Frequency  Domain  Computations 


When  evaluating  roll  motions  in  the  frequency  domain,  an  iterative  procedure  must 
be  used  to  determine  the  correct  roll  amplitude  for  roll  damping  computations.  Con¬ 
vergence  is  obtained  when  the  following  is  satisfied: 


*  /  ~dampinq\  * dampinq 

74(74  )  -  74 


<  eto\f,A) 


(189) 


where  f)4{f)fmipm°)  is  the  computed  roll  amplitude  based  on  roll  damping  for  roll 
amplitude  and  e:/o,(//4  )  is  the  tolerance  on  roll  damping  amplitude.  To 
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ensure  repeatability  of  calculations,  a  small  tolerance  should  be  used  for  roll  ampli¬ 
tude,  and  a  value  of  0.0001  degrees  has  been  selected. 

As  an  illustrative  example,  Figure  14  shows  the  variation  of  computed  roll  ampli¬ 
tude  with  roll  amplitude  used  for  determining  nonlinear  roll  damping.  The  present 
case  has  significant  roll  resonance,  causing  predicted  roll  amplitude  to  be  very  de¬ 
pendent  on  the  roll  amplitude  used  for  computing  damping.  Figure  14  also  shows 
the  line  7)4(174  )  =  r/4  F  .  The  correct  solution  is  at  the  intersection  01  the 

two  lines. 


Figure  14:  Computed  Roll  Amplitude  Versus  Roll  Amplitude  for  Computing 
Damping,  Frigate  in  Regular  Waves 

Iterative  procedures  for  roll  amplitude  often  employ  the  simple  method  of  using  the 
computed  roll  amplitude  from  the  previous  iteration  for  computing  roll  damping. 
This  approach  neglects  the  variation  of  7)4(7] flinptn9)  with  174  when  determining  the 
next  value  for  fj^amPins  _  Much  faster  convergence  can  be  achieved  by  considering 
the  variation  of  1)4(174  )  with  174  ,  applying  the  following  increment  to 

~  dampinq  •  •  .• 

7)4  F  in  a  given  iteration: 
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(190) 


The  term  dfj4(i)fampin9) / di)famptn9  is  evaluated  using  a  secant  method  based  on 
linear  interpolation  of  computed  values  of  fj4(fjfmiinn9)  from  previous  iterations.  If 
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no  values  are  available  from  previous  iterations,  then  /  dr)fimpin9  is 

assumed  to  be  zero. 

Faster  convergence  is  achieved  if  the  first  iteration  is  based  on  a  reasonable  initial 
guess  of  roll  amplitude.  In  long  waves,  the  ship  roll  amplitude  is  given  by: 

fj 4  =  ki  a  |sin/3|  (191) 

where  kj  is  the  wavenumber  of  incident  waves.  In  short  waves,  the  roll  amplitude 
goes  to  zero.  The  following  equation  has  been  developed  based  on  computed  roll 
motions  for  the  initial  guess  of  roll  amplitude: 

574  ~  kj  a  | sin/? |  exp  (—0.08  Lpp  kj)  (192) 


The  above  procedure  has  been  found  to  be  very  stable  and  efficient,  with  conver¬ 
gence  usually  achieved  within  8  iterations.  For  frequency  domain  computations 
for  random  seas,  the  roll  amplitude  for  roll  damping  computations  is  taken  as  1.25 
times  the  RMS  roll. 

8  Validation  of  Roll  Damping  and  Motion 
Predictions  with  the  Haslar  Steered 
Warship  Model 


Experiments  for  the  steered  warship  model  of  Lloyd  and  Crossland  [29]  have  been 
selected  for  validating  roll  damping  and  motion  predictions  from  the  ShipMo3D 
library.  These  experiments  were  conducted  in  the  Manoeuvring  Tank  at  Admi¬ 
ralty  Research  Establishment  Haslar  (now  part  of  the  privatized  QinetiQ).  Lloyd 
and  Crossland’s  experiments  were  selected  for  the  present  validation  because  they 
include  comprehensive  measurements  of  both  ship  motions  and  rudder  motions, 
which  were  driven  by  an  autopilot  with  documented  properties.  The  results  reported 
here  include  conditions  that  were  not  documented  in  Reference  29,  for  which  the 
data  were  provided  directly  to  DRDC  Atlantic. 

Table  2  gives  the  particulars  for  the  steered  warship  model,  and  Figure  15  shows  the 
body  plan.  Tables  3  to  6  give  dimensions  for  the  bilge  keels,  propeller  shaft  brack¬ 
ets,  stabilizer  fins,  and  rudders.  Some  of  the  longitudinal  locations  of  appendages 
in  Tables  3  to  6  have  corrections  from  those  originally  presented  in  Reference  29. 
The  steered  warship  model  included  an  autopilot  with  control  parameters  as  given 
in  Table  7. 
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Table  2:  Main  Particulars  for  Haslar  Steered  Warship  Model 


Length,  Lpp 

5580  mm 

Beam,  B 

671  mm 

Midships  draft,  Trnul 

196  mm 

Trim  by  stem,  ts 

0  mm 

Displacement,  A 

345  kg 

Vertical  centre  of  gravity,  KG 

276  mm 

Metacentric  height,  GM fiuid 

76  mm 

Roll  radius  of  gyration,  rxx 

257  mm 

Pitch  radius  of  gyration,  ryy 

1293  mm 

Yaw  radius  of  gyration,  rzz 

1265  mm 

Figure  15:  Body  Plan  for  Haslar  Steered  Warship  Model 
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Table  3:  Bilge  Keel  Dimensions  for  Haslar  Steered  Warship  Model 


Station  (20  at  AP) 

Span  (mm) 

Root  lateral  offset  (mm) 

Root  above  baseline  (mm) 
Dihedral  angle  (deg,  port  side) 


6.79 

7 

8 

9 

9.69 

68 

68 

68 

68 

68 

239 

239 

243 

248 

248 

107 

107 

82 

64 

64 

-45 

-45 

-45 

-45 

-45 

Table  4:  Propeller  Shaft  Bracket  Dimensions  for  Haslar  Steered  Warship  Model 


Aft  inner 

Aft  outer 

Fore  inner 

Fore  outer 

bracket 

bracket 

bracket 

bracket 

Station  (20  at  AP) 

18.40 

18.40 

17.97 

17.97 

Span  (mm) 

165 

127 

81 

68 

Root  chord  (mm) 

30 

30 

12 

12 

Tip  chord  (mm) 

30 

30 

12 

12 

Root  lateral  offset  (mm) 

21 

132 

72 

157 

Root  above  baseline  (mm) 

148 

163 

103 

125 

Dihedral  angle  (deg,  port  side) 

-47.5 

-99.5 

-42 

-90 

Table  5:  Stabilizer  Fin  Dimensions  for  Haslar  Steered  Warship  Model 


Station  (20  at  AP) 

10.38 

Span  (mm) 

83 

Root  chord  (mm) 

174 

Tip  chord  (mm) 

92 

Root  lateral  offset  (mm) 

273 

Root  above  baseline  (mm) 

75 

Dihedral  angle  (deg,  port  side) 

-46 
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Table  6:  Rudder  Dimensions  for  Haslar  Steered  Warship  Model 


Station  (20  at  AP) 

19.46 

Span  (mm) 

153 

Root  chord  (mm) 

124 

Tip  chord  (mm) 

88 

Root  lateral  offset  (mm) 

79 

Root  above  baseline  (mm) 

178 

Dihedral  angle  (deg,  port  side) 

-83 

Table  7:  Rudder  Control  Properties  for  Haslar  Steered  Warship  Model 


Maximum  deflection  angle 

35  deg 

Maximum  deflection  rate  S'".^er 

35  deg/s 

Deflection  natural  frequency  us 

25.8  rad/s 

Deflection  damping  ratio  Q 

0.85 

Yaw  displacement  gain  k$6 

3.8 

Yaw  velocity  gain  kv5 6 

1.7  s 

Yaw  acceleration  gain  kg6 

0.0  s2 
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For  experimental  sway  motions,  Reference  29  describes  a  procedure  for  conversion 
from  body  axes  to  earth  fixed  axes.  The  present  report  uses  experimental  sway  mo¬ 
tions  without  this  correction  because  it  is  doubtful  that  the  correction  is  applicable 
for  sway  motions  measured  using  an  accelerometer,  as  was  done  during  the  model 
tests. 


8.1  Roll  Damping  in  Calm  Water 

The  Halsar  steered  warship  experiments  included  roll  decay  tests  in  calm  water. 
ShipMo3D  predictions  of  the  experimental  conditions  have  been  made  using  the 
following  methods  for  bilge  keel  drag  coefficients: 

•  Ikeda’s  method  (Equation  (125)), 

•  the  simplified  Kato  method  based  on  roll  amplitude  (Equation  (129)). 

Figures  16  to  23  indicate  that  roll  damping  tends  to  be  underpredicted  when  using 
Ikeda’s  method,  and  overpredicted  when  using  the  simplified  Kato  method.  The 
presented  roll  damping  values  are  normalized  by  the  critical  roll  damping  value 
£>44.  When  considering  which  of  the  two  methods  is  more  accurate  for  predict¬ 
ing  bilge  keel  drag  coefficients,  it  should  be  noted  that  there  is  greater  uncertainty 
with  predicted  roll  forces  as  ship  speed  increases,  in  part  due  to  violation  of  the 
assumption  of  low  speed  when  computing  hull  radiation  forces.  Consequently,  roll 
damping  decay  tests  at  lower  speeds  are  likely  more  appropriate  for  determining 
which  method  provides  more  accurate  predictions  of  bilge  keel  drag  coefficients. 
Figures  16  to  20,  which  include  Froude  numbers  up  to  0.18,  suggest  that  Ikeda’s 
method  is  more  accurate  than  the  simplified  Kato  method  for  predicting  bilge  keel 
drag  coefficients. 
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Roll  Amplitude  774 


Figure  16:  Roll  Damping  Steered  Warship,  Froude  Number  =  0.00 


Roll  Amplitude  774 


Figure  17:  Roll  Damping  Steered  Warship,  Froude  Number  =  0.06 
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Roll  Amplitude  774 


Figure  18:  Roll  Damping  Steered  Warship,  Froude  Number  =  0.09 


Roll  Amplitude  774 


Figure  19:  Roll  Damping  Steered  Warship,  Froude  Number  =0.15 
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Figure  21:  Roll  Damping  Steered  Warship,  Froude  Number  =  0.25 
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Roll  Amplitude  774 


Figure  22:  Roll  Damping  Steered  Warship,  Froude  Number  =0.31 


Roll  Amplitude  774 


Figure  23:  Roll  Damping  Steered  Warship,  Froude  Number  =  0.38 
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8.2  Roll  Motions  in  Beam  Seas  Using  Different 
Methods  for  Bilge  Keel  Drag  Coefficients 

Roll  motions  have  been  computed  in  beam  seas  using  the  following  3  methods  for 
predicting  bilge  keel  drag  coefficients: 

•  Ikeda’s  method.  Equation  (125), 

•  the  simplified  Kato  method  based  on  roll  velocity.  Equation  (127). 

•  the  simplified  Kato  method  based  on  roll  amplitude,  Equation  (129). 

A  primary  motivation  for  the  comparison  was  to  determine  if  bilge  keel  drag  co¬ 
efficients  depend  on  roll  amplitude  as  suggested  by  Ikeda  et  al.  [19],  or  on  roll 
velocity  as  suggested  by  Kato  [14].  Figures  24,  25  and  26  show  roll  motions  in 
beam  waves  for  Froude  numbers  of  0.18,  0.28,  and  0.36.  Figure  24  indicates  that 
Ikeda’s  method  provides  the  best  roll  motion  predictions,  while  the  two  approaches 
based  on  Kato  lead  to  underprediction  of  roll  motions.  The  interpretation  of  results 
at  higher  speeds  is  more  difficult,  in  part  due  to  the  greater  difficulty  in  predicting 
hull  radiation  and  diffraction  forces  as  ship  speed  increases. 


Figure  24:  Roll  RAOs  for  Steered  Warship  with  Different  Bilge  Keel  Drag 
Coefficients,  Beam  Seas,  Froude  Number  0. 18 
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Figure  25:  Roll  RAOs  for  Steered  Warship  with  Different  Bilge  Keel  Drag 
Coefficients,  Beam  Seas,  Froude  Number  0.28 


Figure  26:  Roll  RAOs  for  Steered  Warship  with  Different  Bilge  Keel  Drag 
Coefficients,  Beam  Seas,  Froude  Number  0.36 
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8.3  Motion  Predictions  in  All  Experimental 
Conditions 

Figures  27  to  43  show  comparisons  of  predictions  and  experiments  for  all  avail¬ 
able  experimental  conditions.  The  numerical  predictions  include  ShipMo3D  time 
domain,  ShipMo3D  frequency  domain,  and  the  strip  theory  frequency  domain  pro¬ 
gram  SHIPM07  [8].  The  ShipMo3D  predictions  use  Ikeda’s  method  for  bilge  keel 
drag  coefficients.  SHIPM07  uses  Kato’s  method  for  predicting  bilge  keel  drag 
coefficients. 

The  numerical  predictions  from  both  the  ShipMo3D  library  and  from  SHIPM07 
give  generally  good  agreement  with  the  experiments.  The  excellent  agreement  be¬ 
tween  frequency  domain  and  time  domain  results  from  ShipMo3D  indicate  consis¬ 
tency  between  the  frequency  domain  and  time  domain  simulations.  The  similarity 
of  predictions  from  ShipMo3D  and  SHIPM07  is  likely  due  to  slenderness  of  the 
hull  geometry.  As  expected,  the  numerical  predictions  are  better  for  heave  and  pitch 
than  for  sway,  roll,  and  yaw. 

Agreement  between  predictions  and  experiments  generally  decreases  as  ship  speed 
increases,  likely  due  to  the  assumption  of  low  forward  ship  speed  in  the  radia¬ 
tion  and  diffraction  computations.  The  absence  of  hull  lift  forces  in  the  present 
ShipMo3D  predictions  will  also  affect  accuracy  at  higher  ship  speeds.  The  cor¬ 
respondence  of  trends  for  yaw  and  rudder  predictions  suggests  that  the  rudder  re¬ 
sponse  to  yaw  motions  is  being  correctly  modelled. 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 


Figure  27:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  0  degrees,  Froude 
Number  0.28 


□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 


Figure  28:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  0  degrees,  Froude 
Number  0.37 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 


Figure  29:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  30  degrees, 
Froude  Number  0. 18 


44 


DRDC  Atlantic  TM  2004-227 


Yaw  RAO  \t]q\  /{kja) 


□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 
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Figure  30:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  30  degrees, 
Froude  Number  0.27 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 
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Figure  31:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  30  degrees, 
Froude  Number  0.37 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 
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Figure  32:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  60  degrees, 
Froude  Number  0. 18 
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□  Experiments 
—  3D  time  domain  predictions 
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Figure  33:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  60  degrees, 
Froude  Number  0.27 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 


Figure  34:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  60  degrees, 
Froude  Number  0.36 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 
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Figure  35:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  75  degrees, 
Froude  Number  0. 18 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 


Wave  Frequency  ujf  y/L/g  Wave  Frequency  ujj  y/L/g 

Figure  36:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  75  degrees, 
Froude  Number  0.28 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 
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Figure  37:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  75  degrees, 
Froude  Number  0.36 


52 


DRDC  Atlantic  TM  2004-227 


□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 
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Figure  38:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  90  degrees, 
Froude  Number  0. 18 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 


Figure  39:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  90  degrees, 
Froude  Number  0.28 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 
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Figure  40:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  90  degrees, 
Froude  Number  0.36 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 


Wave  Frequency  uif  \J~L~j~g  Wave  Frequency  ujj  \J  L/g 


Figure  41:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  120  degrees, 
Froude  Number  0.27 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 


Figure  42:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  150  degrees, 
Froude  Number  0.26 
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□  Experiments 
—  3D  time  domain  predictions 

3D  frequency  domain  predictions 
Strip  theory  predictions 
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Figure  43:  RAOs  for  Steered  Warship,  Stern  Quartering  Seas  at  180  degrees, 
Froude  Number  0.26 
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9  Recommendations  for  Future  Work 


The  ShipMo3D  library  is  being  extended  to  model  the  motions  of  a  freely  maneu¬ 
vering  ship  in  both  calm  water  and  in  waves.  Hull  lift  forces  are  very  important 
for  maneuvering  ships  due  to  the  high  angles  of  attack  that  can  occur;  thus,  hull  lift 
forces  should  be  added  to  the  ShipMo3D  library.  Hull  lift  forces  will  also  influence 
oscillatory  motions  in  waves. 

In  the  longer  term,  the  ShipMo3D  library  should  be  modified  such  that  radiation  and 
diffraction  computations  can  model  the  influence  of  hull  forward  speed  on  boundary 
conditions,  thus  leading  to  more  accurate  hull  force  predictions.  This  requirement 
could  be  satisfied  using  a  forward  speed  Green  function  in  either  the  time  domain 
or  frequency  domain. 

10  Conclusions 


The  ShipMo3D  library  now  includes  appendage  and  viscous  forces  for  predicting 
ship  motions  in  waves.  Modelled  appendages  can  include  rudders,  bilge  keels, 
skegs,  and  propeller  shaft  brackets.  Rudder  motions  can  be  modelled  using  a  linear 
autopilot. 

The  ShipMo3D  library  can  include  appendages  in  hull  radiation  and  diffraction 
computations  using  dipole  panels  on  appendage  surfaces.  Computations  for  a  naval 
destroyer  have  shown  that  inclusion  of  appendages  in  hull  radiation  and  diffraction 
computations  has  a  negligible  influence  on  predicted  ship  motions. 

Motion  predictions  have  been  validated  with  model  tests  using  a  steered  warship  in 
regular  waves.  Numerical  predictions  give  generally  good  agreement  with  experi¬ 
ments.  Heave  and  pitch  predictions  are  somewhat  better  than  sway,  roll,  and  yaw 
predictions.  Agreement  between  predictions  and  experiments  tends  to  decrease  as 
ship  speed  increases.  The  current  predictions  assume  low  ship  speed  when  evalu¬ 
ating  hull  radiation  and  diffraction  forces.  Furthermore,  hull  lift  forces,  which  are 
proportional  to  ship  speed,  are  neglected  in  the  current  predictions. 

Near-term  work  will  consider  motions  of  freely  maneuvering  ships,  and  will  include 
hull  lift  forces. 
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Symbols  and  Abbreviations 


A  foil 

Alk 

ij 

Abk 

■rip 

A-X 

a 

acd,  bed 

cip 

d  y  ■  az 

B 

TDbk 

*-*ij 

jjbk 

Tjbk—visc 

-D44 

TDcr 

-D44 

ft  foil 

ij 

ft  foil 
P 

o  foil— vise 
T3  hull— eddy 

£*44 

Bx 

[BSS\ 

[Bdij] 

[Bpp] 

r^bk 

'-'ij 

cd 

r< fou 
^ ij 
rfhull 
^  eddy 
(jhuil 

Clift 

(jrudder 

/~i rudder 
° i6 

[Cm\ 

[ft„] 

[C„  J 

c 

Croot 

Ctip 

Fbk 

pbk—lift 


foil  added  mass 

foil  added  mass  for  motion  perpendicular  to  its  surface 
bilge  keel  added  mass 

added  mass  of  bilge  keel  for  motion  perpendicular  to  its  surface 
hull  cross-sectional  area 
wave  amplitude 

coefficients  for  estimating  bilge  keel  drag  coefficients 

flow  acceleration  perpendicular  to  foil 

flow  acceleration  components  in  y  and  z  directions 

ship  beam 

bilge  keel  damping 

damping  of  bilge  keel  for  motion  perpendicular  to  its  surface 
bilge  keel  viscous  roll  damping 
critical  roll  damping 
foil  damping 

foil  damping  coefficient  for  motion  perpendicular  to  its  surface 

foil  viscous  roll  damping 

hull  roll  eddy  damping 

hull  sectional  beam 

rudder  motion  damping  matrix 

rudder-ship  motion  damping  coupling  matrix 

ship  motion  damping  matrix 

ship  motion  stiffness  from  bilge  keel 

drag  coefficient 

ship  motion  stiffness  from  foil 

hull  eddy  damping  coefficient 

hull  skin  friction  coefficient 

lift  coefficient 

force  stiffness  coefficient  for  rudder  deflection 

rudder  deflection  stiffness  coefficient  for  ship  motion  mode  i 

rudder  motion  stiffness  matrix 

rudder-ship  motion  stiffness  coupling  matrix 

ship  motion  stiffness  matrix 

mean  chord  length 

chord  length  at  root 

chord  length  at  tip 

wave  excitation  force  on  bilge  keel  for  mode  j 
lift  force  on  bilge  keel 
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pbk—visc 
p  foil— vise 
rp  foil— vise 
p  hull— eddy 
phull—skin 

{Fr,} 

Gdipole(x,  xd,  ke ) 
G(x,  xd,ke) 
G0(x,xd,ke) 

Go 


Goo  (*G  'pdi  ke 

GM 

9 

h 

Mh ) 

Ji(h) 

KG 

ke 

ki 

l-a 

5J 

La 
KS  j 

hv 
KS  j 

-kpp 

wu 

wy 

KJ 

^dipole 

^ x  •>  W'z 

77-4 

R 

Rxy 

Rbilge 

Ro(h,v) 

Ri 

Ri(h,v) 

Refoil 


trk 

rl 
r.2 
'  l 


viscous  force  on  bilge  keel 

viscous  force  on  foil 

foil  viscous  roll  moment 

hull  eddy  force 

hull  skin  friction  force 

ship  motion  excitation  force  vector 

Green  function  for  dipole 

Green  function  for  source 

zero  frequency  Green  function  for  source 

frequency  dependent  term  of  Green  function  relative  to  zero 

frequency  Green  function 

infinite  frequency  Green  function  for  source 

metacentric  height 

gravitational  acceleration 

non-dimensional  horizontal  distance  from  field  point  to  source 

Bessel  function  of  the  first  kind  of  order  zero 

Bessel  function  of  the  first  kind  of  order  one 

height  of  centre  of  gravity  above  baseline 

encounter  wavenumber 

incident  wavenumber 

rudder  deflection  acceleration  gain  for  ship  motion  mode  j 

rudder  deflection  displacement  gain  for  ship  motion  mode  j 

rudder  deflection  velocity  gain  for  ship  motion  mode  j 

ship  length  between  perpendiculars 

rudder  motion  inertia  matrix 

rudder-ship  motion  inertia  coupling  matrix 

effective  ship  mass  matrix 

normal  vector  for  dipole  surface 

x,  y,  and  z  components  for  normal  pointing  outward  from  hull 
hull  roll  normal 

distance  from  field  point  to  source 
horizontal  distance  from  field  point  to  source 
hull  bilge  radius 

term  for  evaluation  of  frequency  dependent  Green  function 
distance  from  field  point  to  source  image 
term  for  evaluation  of  frequency  dependent  Green  function 
foil  Reynolds  number 

roll  radius  from  ship  centre  of  gravity  to  bilge  keel 
effective  roll  radius  for  lift  forces 
effective  square  of  roll  radius  squared  for  lift  forces 
effective  cube  of  roll  moment  arm  for  viscous  forces 
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ryy 

f  zz 

S 

3 dipole 
S 

Smax 

Tmid 

Tx 

t-S 

U 

Up 

'LL yi  'LL z 

V 
Vp 

v 

x,y,z 
x,  y,z 
x 

3’aft,  i  Xfore 

Xd 

Xpp 
rjQ  f  oil 

Xs 

a 

P 

r 

A  Ubk 
AUfoil 

^  ^ damping 
^ rudder 

sb 

Xrudder 

°C 

etol(rj 4) 

Cs 

Vj 

Vj 

n 4 

V4 

V 

/X 

p 


roll  radius  of  gyration 
pitch  radius  of  gyration 
yaw  radius  of  gyration 
foil  area 
dipole  surface 
span 

maximum  span 

midships  draft 

hull  sectional  draft 

trim  by  stern 

ship  forward  speed 

flow  velocity  perpendicular  to  foil 

flow  velocity  components  in  y  and  z  directions 

flow  velocity 

flow  velocity  perpendicular  to  foil 

non-dimensional  vertical  distance  from  field  point  to  source  image 
coordinates  in  translating  earth  axes 
coordinates  of  foil  centroid 
location  in  space 

aft  and  forward  x  limits  for  long  low-aspect  ratio  foil 
dipole  location 

x  coordinate  of  ship  forward  perpendicular 

x  coordinate  of  foil 

source  location 

foil  angle  of  attack 

wave  direction  relative  to  ship 

foil  dihedral  angle 

local  flow  velocity  correction  at  bilge  keel 

local  flow  velocity  correction  at  foil 

increment  in  guess  for  roll  damping  amplitude 

rudder  deflection  angle 

hull  boundary  layer  thickness 

command  rudder  deflection  angle 

tolerance  for  roll  damping  amplitude 

damping  ratio  for  rudder  control  system 

ship  motion  displacement  for  mode  j 

ship  motion  velocity  for  mode  j 

roll  amplitude 

roll  velocity  amplitude 

water  kinematic  viscosity 

dipole  strength 

water  density 
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local  roll  velocity  ratio 
(f>  velocity  potential 

luc  wave  encounter  frequency 

c os  natural  frequency  for  rudder  control  system 

A  ship  mass  displacement 
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